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Abstract 

We study an initial boundary value problem on a ball for the heat-conductive 
system of compressible Navier-Stokes-Fourier equations, in particular, a criterion 
of breakdown of the classical solution. For smooth initial data away from vacuum, 
it is proved that the classical solution which is spherically symmetric loses its 
regularity in a finite time if and only if the density concentrates or vanishes or the 
velocity becomes unbounded around the center. One possible situation is that a 
vacuum ball appears around the center and the density may concentrate on the 
boundary of the vacuum ball simultaneously. 

I 

1 Introduction and main results 

We are concerned with the heat-conductive system of compressible Navier-Stokes- 
Fourier equations which reads as 

'pt + V-(j>U) = 0, 

< (pU) t -I- V • (pU <g)U) + VP = pAU + (p + A)V(V • U), (1.1) 

k cy ({ P 0)t + V • (pdU)) + PV-U = k A6 + T[V[/], 

where ^ 

T[W] = 2p(V(U)) 2 + \(V -U) 2 , V(JJ) = VU+ 2 W (1.2) 

and t > 0,x € O C M. n (N = 2,3), p = p(t,x),U = U(t,x ) and 0 = 9{t,x) are the 
density, fluid velocity and temperature respectively. P = P(p, 9) is the pressure given 
by a state equation 

P(p) = RpQ. (1.3) 

The shear viscosity p , the bulk one A and heat conductivity n are constants satisfying 
the physical hypothesis 

N 

p,K>0, P + —X>0. (1.4) 

The domain Q is a bounded ball with a radius b, namely, 

H = Bb = {x G M. n -, |x| < b < oo}. (1.5) 

We study an initial boundary value problem for (1.1) with the initial condition 

iP,U,9)(0,x) = (p 0 ,U 0 ,e 0 )(x), xeQ, (1.6) 
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and the boundary condition 


U = 0, — = 0 x £ <90, n is outnormal vector. (1-7) 

on 

We are looking for the smooth spherically symmetric solution (p, U) of the problem 
(1.1), (1.5),(1.6) which enjoys the form 

p(t,x) = p(t,\x\), U(t,x) = u(t,\x\)-^-r, 9(t,x) = 9(t,\x\). (1.8) 

\x\ 

Then, for the initial data to be consistent with the form (1.7), we assume the initial 
data (po> Uo) also takes the form 

Po = Po(\x\), U 0 = u 0 {\x\)-^~, 9 q = 9q{\x\). (1.9) 

\x\ 

In this paper, we further assume the initial density is uniformly positive, that is, 

po = p 0 (|z|) > P > 0, xGfl (1.10) 


for a positive constant p. 

Here, it is noted that since the assumption (1.7) implies 

U(t,x)+ U(t,—x) = 0, iGll. (1.11) 

We necessarily have U(t,0) = 0 (also Uo{0) = 0) as long as classical solutions are 
concerned. 

There are many results about the existence of local and global strong solutions in time 
of the isentropic system of compressible Navier-Stokes equations when the initial density 
is uniformly positive (refer to [Hl9l- llllll9I[2Tll24ll25j and their generalization [I5HI71I23] 
to the full system including the conservation law of energy). On the other hand, for 
the initial density allowing vacuum, the local well-posedness of strong solutions of the 
isentropic and heat-conductive system was established by Kim mm- For strong 
solutions with spatial symmetries, the authors in [14] proved the global existence of 
radially symmetric strong solutions of the isentropic system in an annular domain, 
even allowing vacuum initially. 

However, it still remains open whether there exist global strong solutions which are 
spherically symmetric in a ball. The main difficulties lie on the lack of estimates of 
the density and velocity near the center. In the case vacuum appears, it is worth 
noting that Xin (26) established a blow-up result which shows that if the initial density 
has a compact support, then any smooth solution to the Cauchy problem of the full 
system of compressible Navier-Stokes equations without heat conduction blows up in a 
finite time, see more generalization in m, m .The same blowup phenomenon occurs 
also for the isentropic system. Indeed, Zhang-Fang ( [28] .Theorem 1.8) showed that 
if ( p,U ) £ C' 1 ([0, T]; H k ) (k > 3) is a spherically symmetric solution to the Cauchy 
problem with the compact supported initial density, then the upper limit of T must be 
finite. To deal with large discontinuous data, Hoff [3] established global weak solutions 
of the symmetric compressible heat-conductive flows. On the other hand, it’s unclear 
whether the strong (classical) solutions lose their regularity in a finite time when the 
initial density is uniformly away from vacuum. Therefore, it is important to study the 
mechanism of possible blowup of smooth solutions, which is a main issue in this paper. 


2 


In the spherical coordinates, the original system m under the assumption (1.8) 
takes the form 


where 


Pt + ( pu )$ = 0, 

(pu)t + (pu 2 )t +P r = 

c v ((p9) t + {pu0)() + Pu£ = n9 r £ + v(u{) 2 - —(r N ~ 2 u 2 ) r , 


d d N -1 

v — 2p + X, — — — H . 

o^ or r 


Without lossing of generality, we assume cy = 1 and N = 3. 
Now, we consider the following Lagrangian transformation: 


pr 

t = t. h= p(t,s)s 2 ds , rj = (pr 2 ) 

Jo 

Then, it follows from (1.12) that 

, u 

h t = -, n = u, r h = p , 

V 

and the system (11.121) can be further reduced to 

( r 2 r])t = ( r 2 u ) h , (<S=> = u h ), 

2 


2 \ —1 


u t = r ( —R—^— + v 

rjr* rjrj 


9 


u h , 2 

-1 u 

r h r 


u t + Rr~ -TT- = vr 

\ r vJ h 

r 2 u) h 


9 t = -R9- 


= -R9 


ry* 2 

[r 2 u)h 


{r 2 u)i 

r 2 p J h 

u h , 2it\ 2 


+ ur 2 p —- H-— 4 p{ru 2 )h + k 


rh r J 
Uh , 2 


r 2 9 h 


r h J h 


+ A r-p -1- -u + 2pr rj H- T + k 


u\ 2 u 2 


r 2 r 2 
'h ' 


Jh r 

The initial boundary value problem for system (11.161) 

(u,r],9)(0,h) = (u 0 ,rio,9o), (r/ 0 > 0,9 0 > 0), 

u(t, 0) = u(t, 1) = 0, 9 h (t, 0) = 9 h (t, 1) = 0, 

where t > 0, h E [0,1] and 


r 2 9 h 
r h J h 


1 = 


["pot 

Jo 


r)r 2 dr = / p(t,r)r 2 dr , 


according to the conservation of mass. Note that 

r{t , 0) = 0, r(t , 1) = b. 

We denote Eq the initial energy 

Eo = J |y + R (jIpq - log t-q^o - 1) + (0 O - log 9q - 1)| dh. 

Our main result is stated as follows. 


( 1 . 12 ) 


(1.13) 


(1.14) 


(1.15) 


(1.16) 


(1.17) 


(1.18) 


(1.19) 


( 1 . 20 ) 
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Theorem 1.1 Assume that the initial data (pq,Uq,9 0 ) satisfy (1.8),(1.9),(1.10) and 

(po,U 0 ,6 0 ) GH 3 (n). (1.21) 

Let (yO, U, 9) be a classical spherically symmetric solution to the initial boundary value 
problem (1.1),(1.5),(1.6),(1.7) in [0,T] x O, and T* be the upper limit ofT, that is, the 
maximal time of existence of the classical solution. Then, if T* < oo, it holds that 

lim sup ( p(t, |x|) + -(t, |x|) + \U\(t, |x|) ) = oo. (1-22) 

(t,|ai|)—>(T*, 0 ) V P / 

Remark 1.1 The local existence of smooth solution with initial data as in Theorem 
o is classical and can be found, for example, in m and references therein. So the 
maximal time T* is well defined. 

Remark 1.2 There are several results on the blowup criterion for strong and classical 
solutions to the isentropic and heat-conductive system HI.ID (refer to [41471 1221129] and 
references therein). Especially, the authors in [?| established the following Serrin-type 
blowup criterion: 


lim sup (||p||l°°( 0 ,T;L°°) + \\U\\ L r(o, T;L«)) = oo, (1-23) 

T/T* 

for any r £ [2, oo] and s £ (3, oo] satisfying 

- + -<1. (1.24) 

r s 

Theorem 1 1.1 1 asserts that the formation of singularity is only due to the concentration 
or cavitation of the density and velocity around the center. More precisely, the density 
anywhere away from the center is bounded up to the maximal time. Also recall that 

U(t, 0) = 0, for all t < T* (1-25) 

as far as classical solution is concerned. It indicates the possible loss of regularity of 
velocity at the center. 

Remark 1.3 Theorem 1.1 may be viewed as an extension of recent work in [8] 
where the authors established blowup criterion for baratropic spherically symmetric 
Navier-Stokes equations. 


2 Proof of Theorem 11.11 

We only prove the case when N = 3 since the case N = 2 is even simpler. Throughout 
of this section, we assume that (p, U, 9) is a classical spherically symmetric solution with 
the form (1.8) to the initial boundary value problem (1.1),(1.5),(1.6),(1.7) in [0, T] x 0, 
and the maximal time T*, the upper limit of T, is finite, and we denote by C generic 
positive constants only depending on the initial data and the maximal time T*. 

By simple calculations, we have the following estimates. 
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Lemma 2.1 It holds that, 


/ r/dh = / r/odh , 
Jo Jo 


J r2 vdh = ^ ( r 2 r/) 0 dh = y, 


iw 2 + 0 ) dh = 


+ 0o ) dh. 


Also, we have the following basic energy equality. 
Lemma 2.2 It holds that 

£{t)+ [ V(r)dr = £( 0), 


where 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


(2.4) 


£(t) = j |y + R (r 2 rj — logr 2 rj — l) + (6 — log# — l)j> dh, 


m = 


Xr 2 r] f uh 2 


n i I u 1+2 „ 

9 \r h r ) 0 


fir 2 r) I 2 u 


— + 

r 2 ^ ,„2 

'h ' 


(+ + o aO 


2 r 2 r/ / Uh 2 


# W 


+ -it + 


+ re 

2 _2— /-. \ 2 


r20|) 

Di # 2 J 


dh 


+ re^ 
0 V " r 'h rj r h 9 2 


fir r) Uh u 


dh, 


or equivalently, 


/ 

Jo 


pSdx + 



o Jn 


where 


1 |v0| 2 ' 

-\h[VCJ] + re ) dxds = / p 0 S 0 dx, 


1 


5 = i^Oo" 1 ) + 4>(0) + -|C/| 2 , 4*(s) = s - logs - 1. 


(2.5) 


( 2 . 6 ) 


(2.7) 


In order to prove Theorem 1.1, we can show the following stronger characterization 
of blowup criterion, that is, 

lim sup (p{t,h) + ~{t,h)+ \\U\\ L 2 u T *. L ^ {Bh)) ) = oo. (2.8) 

(t,h)—>(T* — ,0) \ P J 


We argue by contradiction. For original system (1.1), assume there exists a small 
r\,e > 0 and a constant C, such that 

p(t,r) + ~{t,r) + \U\(t,r) < C, for (t, r) € (T* - e, T*) x [0, n]. (2.9) 

P 

Through the Lagrangian transformation (11.141) , one immediately conclude that for sys¬ 
tem (|1.16p . there exists a small constant hi > 0, such that 

p(t,h) + i(t,h) + \U\(t,h) < C, for (t, h) € (T* — e, T*) x [0, h\], (2.10) 

Denote 

ho = -h\. ( 2 - 11 ) 

Recall blowup criterion (11.231) by taking r = oo, s = 2, it amounts to prove the following 
Proposition. 
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Proposition 2.3 For ho = \h\, there exists a constant C depending on ho such that 
p(T,h) + n| L 2 ( 0 ,T;L~(Be)) < C(ho) for ( T,h ) G (T* — e,T*) x [h 0 , 1], (2.12) 


To do that, we prepare the next lemma which gives a relationship between r and y. 


Lemma 2.4 There exists a positive constant C independent of T and two strict in¬ 
creasing function 0j : [0,1] —> [0, oo) such that 


r(t,h) > C0i(/i) 

and 

b 3 — r(t, h) 3 > CQ 2 (h), 


for all 0 < t < T*. 

Proof. For s > 0, set 

G(s) = slogs — s + 1. 

Obviously, G is a convex function in (0, oo). By Jensen’s inequality, one has 


G 


J Br P dx \ / !b t G(p)dx 


I B r 


< 


I B r 


G 


/CVi\ ^ Ci f n G(p)dx 
\ r 3 J ~ r 3 


(2.13) 

(2.14) 


(2.15) 


(2.16) 


Consequently, the uniform estimates for r(t,h) follows immediately from Entropy in¬ 
equality (|2.6p and (12.16[) . That is, given h > 0 there is a strict increasing function 
0i (h) such that 


r(t,h ) > 0i (h) 0i(O) = 0 and ©i (h) > 0 for h > 0. (2-17) 

By a similar step, one can conclude (|2.14l) by the following Jesen’s inequality 


G 


' J B-P dx \ lB- G (p) dx 


IBS: 




IBS 


(2.18) 


We are now in a position to establish the pointwise estimates of the density away 
from the center. To do that, we first write the density in the following form. One may 
refer to [T8J for a similar representation in one-dimensional case. 

Lemma 2.5 

^ h) = BxmtJ) ( (r2,?>o< ' i)+ 1 ? B(t ' h)y(T ' h)0(T ■ h)dT ) - <2J£l) 

Here 

B(t, h) = exp -{ [ k ^df - [ k ^di - ( I' r 2 V dh 
v l Jho r o Jh 0 r \Jh 0 

x f f (r 2 ry) 0 / ^d^dh — f (r 2 r/) f d£dh + v f ( r 2 rj)dh — v f ( r 2 r])odh 
\J ho Jho r 0 J ho J ho x J ho J ho 

( 2 . 20 ) 
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and 




-1 



u 2 + R9 + r 2 r] 


ho 


r h 2 u 2 

An r 3 


df) dh + 


ho 


(r 2 u) h dh 



a(s, ho)ds 


where 


{r 2 rj) t 


a(t,h) =-R—~-+ v- „ 

ry* Zj /yj 


Proof. In view of (11.161) . it holds 


1 / 9 

~^u t = -R-z- + v 

l ry* ryj 


( r 2 il) t \ A 


= ( <r(t,h)) h . 


r V J h 

Thus, for h > ho > 0, integrating (12.23|) over (ho, h ), we deduce that 


ho 


rh 2 u 2 


~ 2 d £) + / ~^rdf = a(t,h) - a(t,h 0 ). 


' ho 


Multiplying 4-2. Q n both sides of (|2.24l) yields 


2 

r rj 



u 


ho 


df) + 


r h 2 u 2 


t J h 0 


R. 


df + aft, h 0 ) > = - 6 + (r 2 rj) t , 


which is 


Denoted by 


Aft, h) = - 



ho 


In view of (|2.26l) and (12.271) . one has 


fr 2 rj)(t, h) = exp 
On the other hand, recall (12.23[) ~ (|2.24D 

rh 


a(t,ho) = - 


ho 


df) ~ 


f 

J ho 


h 2 u 2 


6 U 2 . 


(r 2 ri)t 

T d£-R-z- + v^±. 

3 r Zrjj r Ar^ 


Multiplying r 2 ?/ on both sides of (|2.29l) . the hrst term becomes 


-( r 2 V ) / -odf) = - 


ho 


(r 2 r,) I" ^d(X + (r 2 r,) t P ^d£ 

Jho r J t ho r 2 

eh 


U 


1 ho 


= ~U r v) ^ + (»'«)/ -o d f> -(r u)-o ■ 


1 ho 



( 2 . 21 ) 

( 2 . 22 ) 


(2.23) 


(2.24) 


(2.25) 


(■ r 2 r))t --{( [ ~ 2 ( h'] + [ ‘~pdf + cr(t,h 0 )| fr 2 rf) = —0. (2.26) 

v l \Jho r Jt Jho r J v 


—df) + [ —pdf + aft, ho)l . (2.27) 

r ) t Jh 0 r J 


J Adr^j ■ |(r 2 ?y)o(h) + J exp (^J Ads^j ■ — . (2.28) 


(2.29) 


(2.30) 
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Integrating (12.291) on [/i 0 ,1] x [0, t\ and taking into account the boundary conditions 
(11.171) yields 


n rV( r i ho)dhd 

o 


t r l 


T = — 


'0 J ho 

rt r l 


(r 2 ri) / —d£ \ dhdr + 


' ho 


t r i 


( r 2 u ) 




I£{“ 2+ ™ +(r2,,) £^ 


0 J ho 
d£ > dhdr 


' ho 


' ho 

rt r l 


+ v 


[ [ (r 2 n), 
J 0 J ho 


dhdr. 


(2.31) 


Therefore, 

rt r l 


r 2 r]a(T,h 0 )dhdT = / (r 2 r ?) 0 • 


0 J h 0 


'ho 
rt rl 


/ ^ 

'h 0 r 0 


di ) dh ~H 



I ho 

|ri 2 + R0 + ( r 2 r]) J — ^-df, | dhdr 

f ( r 2 rj)dh — v f ( r 2r q)odh. 

Jho Jho 


+ V 


(2.32) 


In view of boundary condition and (II.I 6 D 1 ,the l,h.s of (12.321) can be written as 


f(f 

JO \Jh 0 
»1 


r 2 gdh 


a(s,ho)ds) dr 


'ho 

f 

'h 0 


r 2 r]dh I ( / cr(r, ho)dr I — 


t / rl 


r 2 rjdh 


/ a(r,h 0 )di 

Jo 


tO \Jho 
rt / rho 


r + 


( r 2 rj) t dh 
rtio 

J {r 2 u) h dh 


a(s 1 ho)ds ) dr (2.33) 


cr(s, ho)ds I dr 


Collecting (I2.28D - (I2.33I) . the proof of Lemma 12.51 is completed. 

We immediately have the following corollary 

Corollary 2.6 Given 0 < ho < 1. for h$ < h < 1, there exists a constant C depending 
on ho such that 

C~\ho) < B(t , h),y(t, h) < C(h 0 ), (2.34) 


and 


C 1 (ho)exp{C' 1 (ho)(t — t)} < ^ 77 —yy < C(ho) exp{C'(h, 0 )(t — r)}, 0 < r < t. 


y(r,h) 


Proof. First, 


/ r 2 r]dh = / —dh = / r z dr = 
J ho Jho P Jr\ 


b 3 — r(ho ) 3 


'r{h 0 ) 


In view of (12.14[) . one gets 


0 < C0 2 (/i o ) < f r 2 i]dh < b 3 . 
J ho 


(2.35) 

(2.36) 

(2.37) 


dhdr 


















Similarly, 


0 < 


io rS Jr (, 


pu 2 r 2 


dr < Cr{h 0 )~ 6 < C0i(/i o ) _ • 
1/2 


lh 0 T “ Jr{h 0 ) r 

JJ^dhl < min r(t, h)~ 2 [ \u\dh < C@i(ho)~ 2 ( [ u 2 dh\ (1 — ho) 1 / 2 

I h 0 r 2 h>h 0 J hQ \J o J 

< ceM- 2 . 

Also, one can verify 


r i r h 9 ?/ 2 
/ 7,2 V / -^rd^dh 
I ho Jh 0 r3 


is bounded. 

Observing that 


f 1 2 [ h 2u 2 

Inr’ln,— 


1 ( /<*.3 rh o„,2 


d£dh = 


/" 2rd : \ r 3 2w 2 . 

:„,^L- d V h -~ Uh 

b 3 f 1 2u 2 ,. [ l 2u 2 „ 

-a£ — / — dh 

Jh, 


3 Ao r 3 


'ho 


is bounded from below and above by a constant C(ho). 
To finish the proof, it suffices to bound 

rt 


<t(t, ho)dr. 


>o 


(2.38) 

(2.39) 

(2.40) 


(2.41) 


(2.42) 


Recall that the r.h.s of (J2.32I) is bounded by some constant C(ho) and (12.331) . one has 

rt pt / rtiQ 


ft ft / fho 

/ a(T,h 0 )dr\ < C(h 0 ) + C(h 0 ) / (| / (r^)^^ 

4o Jo \ Jo 


/ er(s, ho)ds\ ) dr. , 

I o 7 (2-43) 


On the other hand, with the help of (12.101) yields that 

fh 0 _ fho 

Jo 


) fh o 

(r 2 u) h dh = / ( r 3 -) h dh 

Jo r 


rho 


rho 


/ 3 rpudh + / r 6 — dh 

1 0 Jo KrJ h 


u 


<C max |r 2 77| 2 ( / ridh 

0<h<ho V./o 


2 / r^o \ 2 

. 2 , 


u dh 


+i ( 7 (f 

/ 0 rjO \r J h 


fho 

0 
».4 


r 2 p0dh 


< 




1 

2 \ 2 
h 


= C(ho) + C(h 0 ) 


1 

r^° r 2 ^ / Uh u \ 2 \ 2 


/o 


0 \r h r 


(2.44) 


<C'(/t 0 ) + C'(/ l 0 )V(t) 2 , 
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where we used energy inequality, (12.11) , (12.31) and (12.461) in the next Lemma 12.71 

Consequently, the desired bound for <r(r, ho) follows immediately from Gronwall’s 
inequality and (I2.43l) - fl2.44[) . 

Thus finishes the proof of Corollary 12.61 

Hence, substituting (12.201) and (12.211) into (12.191) . we finally arrive at 


Lemma 2.7 


and 


C i (ho) < p(t,h) < C(h 0 ), 0 < h 0 < h, 


C~ l < I Odh < C, 


/ max u 2 (r,h)dr < C(ho). 

I o h&[ho,l] 


(2.45) 

(2.46) 

(2.47) 


Remark 2.1 Theorem 1.1 follows immediately from Lemma 12.71 

Proof. The right-hand side of (12.451) is a direct consequence of the fact that 9 > 0, 
(12.191) and (12.341) . It remains to show the upper bound of r 2 r/ = j y 

Step 1. Multiplying l on both sides of (|1.16)L yield 


(log 0) t = -R(log(r 2 i])) t + (A + |/i)^p 


. 1 r 2 e h 
+ K I a ■ - 

o r h 


dh r 2 0 h 

+ -' 

9 2 r h 


u h , 2 

- 1 --u 

r h r 


+ 


pr 2 r] 

~9~ 



(2.48) 


Integrating (|2.48l) over [0,1] x [0, t] and recall (1.4), (11.171) to obtain 

{ f log Odh — I log 9odh} < —R{ [ log (r 2 r])dh — f log(r 2 r/)od/i}. (2.49) 

Jo Jo Jo Jo 

Applying Jesen’s inequality to (|2.49l) and (2.2) to arrive 


log 


[ Odh > f 
Jo Jo 


> / log Odh 


> f log Oodh — R I log (r 2 r])dh + R I log(r 2 r])odh 

Jo Jo Jo 


-1 r l 


(2.50) 


f { log^o + Rlog(r 2 ri)o}dh, — Rlog f r 2 r)dh 

Jo Jo 


> 

/0 
= C U 


which gives the desired bound for fj Odh. 

Step 2. Applying the mean value theorem to continuous function Oft, h) to get 


f, 0(t,h)dh 

Wt > 0, 3h{t) £ [h 0 ,l], s.t 0(t,h(t )) = -2----. (2-51) 

1 — ho 
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Therefore, for h > ho 


9(t,h)z = 9(t,h(t))z + [ - ——— r d£ 

Jh( 


h(t) 20(i,£) 2 


< C(1 — h 0 ) 


[ 9dh )3 

Jo 


+ 


'ho 

f 1 9 2 

<c{ 1+ max r 2 r](t,h ) / 

1 hG[h 0 ,i] J ho r 2 r]9 2 

<C(/io){l+ max r 2 r)(t, h) ■ V(i)} 
he[h 0 ,l] 




1 

h 


'ho 


49 2 


(2.52) 




Consequently, 


0(t,/t) < C(/io){l + max r 2 7j(t, h) ■ V(t)|, ho < h < 1. 


(2.53) 


S7ep 5. Observing that 
r 2 r/(t, h) = 

< 


B(t,fc)W,>0 + 1 *B(.T,h)y(T,h)0(T,h)dT 


1 


B(t,h)y(t,h) 


{r 2 r])o(h) 


(2.54) 


Hence, 


max r 
fee [ho,i] 


+ C /‘5hl§.ZCfi(l + max rVt,ft)-V(t) 
Jo B(t,h ) T(i,h)\ he[h 0 ,i] 


l rj(t,h) < C + C f exp{—ajf — r)}/l + max r 2 r/(r, h) ■ V(r)ldr 
Jo l he[h 0 ,i] J 

C + C exp{—a(t — r)} max r 2 r)(r, h) • V(r)dr. 

Jo hG[h 0 ,l] 


< 


Write 

£jf) = / exp{—a(i — r)} max ?’ 2 ??(r, /r) • V(r)dr. 

Jo he[h 0 ,l] 

One immediately has 

E t < max r 2 r](t , h) ■ V(t) < (C + CE)V(t ) — aE, 
he[ho,l] 

Tit + (a - CV(i))£7 < CV(f), 

Applying Gronwall’s inequality to (12.571) yields 

E < C exp | — J a — C'V(r)dr| x C J exp | J a — CV{s)ds^V(r)dT. 

The upper bound of r 2 i] follows from ([2.551) and (12.581) . 

Step 4 • It suffices to establish a bound for the velocity. Indeed, 


(2.55) 


(2.56) 


(2.57) 


(2.58) 


(2.59) 
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That is 


max (-) (h) < C(h 0 )V(t) € L 1 (0,T). (2.60) 

/ie[/i 0 ,i] Vr/ 

To conclude, (12.471) is a direct consequence of (12.601) . 

This finishes the proof of Lemma 12.71 
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